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LERWICK ANEMCGRAPH RECORDS 1957-70 AND THE 
OFFSHORE INDUSTRY 


By H. C. SHELLARD 


Summary. Hourly mean winds from Lerwick Observatory for the 14 years 1957-70 are 
analysed to show the frequencies of strong winds of various durations and speeds from different 
directions. Details of some outstanding storms are given and estimates are made of probable 
extreme speeds for various averaging periods. The results, modified to represent conditions 
over the open sea, are used in conjunction with some well-known wave/wind relationships to 
provide estimates of wave-height frequencies for different direction ranges and of probable 
extreme wave heights in the Shetland area. 


Introduction. Increasing mineral exploration in the northern North Sea in 
recent years with the discovery of some important oil fields to the east of the 
Shetland Islands has led to an increased demand from the offshore industry for 
detailed wind analyses for use in the design and operation of offshore structures. 
In order to answer such questions as ‘how often is a wind of 40 kt or more likely 
to blow for 12 hours or more from a given direction?’ or ‘what is the wind speed 
averaged over 6 hours that is likely to be exceeded only once in 50 years?”, it is 
essential to have a long series of continuous wind records in computer-accessible 
form. Such observational material is simply non-existent for the offshore sea 
areas, nor can it possibly be made available for at least another 10 years, by 
which time the need for it will probably have become less urgent. 

However, we already have long-period continuous wind records for Lerwick 
Observatory and the data for the 14 years 1957-70 are available on magnetic 
tape. These data can be analysed to answer the sorts of question mentioned 
above and the results used to provide good estimates of wind conditions over the 
open sea in the general vicinity of the Shetland Islands, together with an indica- 
tion of the wave conditions that they are likely to generate. Attention will be 
mainly confined to wind speeds averaging 25 kt or more and to durations of 
between 1 hour and 24 hours. 


Frequencies of storms of various durations. Table I gives for each of 
the 12 30-degree direction ranges (350-010°, 020-040°, — — — — 320-340°) the 
numbers of spells of various durations during which successive hourly mean 
speeds were 25 kt or more, 30 kt or more and so on. The first part (350-010°) 
shows, for example, that there were six occasions during the 14 years when for 
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spells of g to 11 hours all hourly mean speeds were 30 kt or more from a direction 
in the range 350-010°. It should be noted that the frequencies are non- 
cumulative with respect to duration, i.e. a spell of, say, 12 hours with speeds 
of 30 kt or more is counted only once, in the column headed 12-14 hours, and is 
not also counted as two six-hour spells, four three-hour spells and so on. Hence, 
addition of the number in column 1, three times the number in column Q, six 
times the number in column 3, and so on, will approximately give the total 
number of hours with speeds equal to or greater than the speed concerned. 
However, exact total hours are given in the final column of each table. Inci- 
dentally, if these values are divided by the total numbers of spells given in the 
previous column, average spell-lengths for each speed level and direction are 
obtained and it is of some interest that these tended to be longest for wind 


directions from between 110 and 190° and to decrease in length with increasin 
. S 


speed. 

An important conclusion from Table I is that the most severe storms (mean 
speeds of 55 kt or more throughout) came fron: directions between 200 and 280°. 
However, if attention is confined to storms having durations of three hours or 
more and mean speeds of 40 kt or more, that is to say, to storms which are likely 
to produce significant wave heights of over 5 metres, assuming that fetches are 


not limiting, their distribution with wind direction was as follows: 


Direction Number of Direction Number of 
range storms range storms 
350-010 ‘ 170-—190° 7 
020-040 ‘ 200-220° 6 
050-070° 230-—250° QI 
080-100 260-280° 13 
110-130° 290-310° 3 
140-160 ) 320-340° 2 
Total 64 


Clearly, storms which are of sufficient intensity and duration to produce high 
waves may come from a wide range of directions in the Shetland area. How- 
ever, more than half of them came from between 230 and 280°,while directions 
between 050 and 130° produced none and so seem unlikely to develop very big 
seas. The anemograph at Lerwick has a good open exposure in all directions. 
The only shelter from easterly winds is provided by the small island of Bressay, 
bearing about 030 to 110°, rising to a height of about 225 m at a distance of 
about 53 km. Its effect is not likely to be very great as its terrain is relatively 
smooth and treeless. 

In the original computer tabulations, frequencies were included for every 
duration, at 1-hour intervals, from 1 to 35 hours and for the ranges 36-41, 42-47, 
48-59, 60-71 and 72 hours or more. They were given separately for the periods 
April to September and October to March as well as for the year as a whole. To 
give an indication of the seasonal variation the total frequencies for all directions 
are given in Table II for ‘summer’ (April—Sept.), ‘winter’ (Oct.-March) and for 
the whole year. Although the winter half-year accounts for a large majority of 
the storms, over 70 per cent in most categories, the summer half-year does 
experience an appreciable number with mean speeds of 40 kt or more lasting for 
several hours. However, every one of the 13 storms with mean speeds of 45 kt or 
more throughout, which occurred in the summer half-year took place either in 
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April or in September, there being none at all during the months May to August 
over the 14-year period. This is shown in Table III which gives the monthly 
distributions of the number of hours with mean speeds of 25, 30, 35, 40, 45, 50, 
55 and 60 kt or more, irrespective of direction. It may also be seen that over the 
14 years there were only 14 hours, one per year on the average, having mean 
speeds of 40 kt or more during the four months May to August, compared with 


1 


491 such hours in the other eight months of the year. 


TABLE I—FREQUENCY OF SPELLS OF DURATIONS FROM I HOUR TO 24 HOURS OR 

MORE WITH HOURLY MEAN WIND SPEEDS EQUAL TO OR GREATER THAN THE STATED 

VALUES AND WITH MEAN IND DIRECTIONS WITHIN THE STATED 30-DEGREE 
RANGES AT LERWICK (Ni sERS OF OCCURRENCES IN THE I4 YEARS 1957- 79) 


Duration (hours) 
6-8 g-1I! 12-14 
Numbers of occur 


350-010 


“ 
= OI Nd 
Nn wens Or© 


16 II : I 
6 6 


2 


020-040° 
4 


2 


050-070° 


o80—100° 
9 


110-130° 
I 
I 


140-160° 
3 
3 
I 


170—1g90° 
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TABLE I—continued 
Duration (hours) 
6-8 g-Ir 12-14 15-17 18-23 >24 
Speed Numbers of occurrences 
(kt) > 200-220° 


25 53 14 
30 hy 13 4 
35 3 2 
40 
45 
50 
55 
60 


230—250° 
8 


25 17 
30 ; = 3 I 
33 S! I 
40 
45 
50 


55 
60 


25 
30 
35 
40 
45 
50 
55 
60 


290-310° 
25 : 6 2 
30 
35 
40 
45 
50 

320-340° 

25 3 
30 ; 66 
35 25 
40 s 7 
45 : 3 
50 I 


Details of some outstanding storms. Questions are raised from time to 
time about the profiles or structure of typical wind storms; for example, how 
long do they last, how quickly are their maximum speeds attained, how soon do 
they die away and is there a typical storm profile which could be assumed for 
engineering purposes. 

In Figure 1 the hourly mean speeds and directions in each of the 15 storms in 
which speeds reached or exceeded 50 kt at Lerwick during the 14 years 1957-70 
are plotted against time. All hours in which mean speeds were 30 kt or more are 
included, that is to say, 30 kt is taken as the threshold speed defining the begin- 
ning and end of each storm. This is a quite arbitrary rule, but some such rule is 
necessary as a basis for answering the questions mentioned in the preceding 
paragraph. On this basis the durations of the storms averaged abcut 17 hours, 
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ranging from 6 hours (number 7) to 37 hours (number 6). If the threshold speed 
had been set at only 25 kt the durations would have averaged about 25 hours, 
ranging from 15 hours (number 1) to 54 hours (number 2). It can be seen that if 
the threshold speed had been 40 kt then durations would have ranged from only 
2 hours (number 13) to 21 hours (number 6). In some of the storms the maxi- 
mum speed was reached in a few hours, for example in numbers 7, 10 and 13, 
but in others much more slowly, for example numbers 3 and 6. The rates of 
decline were also very variable. In only one of the 15 storms (number 6) did the 
wind direction back with time; in all the others it veered, sometimes gradually, 
sometimes more sharply. The total direction change, for a 30-knot threshold 
speed, ranged from about 20 degrees to about 100 degrees. 

There is an apparent discrepancy between the 15 storms of Figure 1 and the 
20 spells with hourly mean speeds of 50 kt or more shown by Table II. The 
explanation lies in the fact that some of the storms of Figure 1 include more than 
one of the spells listed in Table II, which represent all occasions of one hour or 
more with speeds of 50 kt or more from within one of the 12 fixed 30-degree 
direction ranges. Thus storm number 3 (Figure 1) includes a one-hour spell 
from 200-220° and a two-hour spell from 230-250° while storm number 6 in- 
cludes one three-hour and one two-hour spell, each from 350-010° but separated 
by an hour with mean speed below 50 kt. 

The storms whose profiles are shown in Figure 1 were selected on the basis of a 
high maximum speed, irrespective of wind direction. Some of these may not 


TABLE II—FREQUENCIES OF SPELLS OF DURATIONS FROM I HOUR TO 24 HOURS OR 

MORE WITH HOURLY MEAN WIND SPEEDS EQUAL TO OR GREATER THAN THE STATED 

VALUES AND FOR ALL DIRECTION RANGES COMBINED AT LERWICK (NUMBERS OF 

OCCURRENCES IN THE I4 YEARS 1957-70 FOR SUMMER (APRIL—SEPTEMBER), 
WINTER (OCTOBER—MARCH) AND YEAR) 


Duration (hours) 

is 3-5 6-8 g-11 12-14 15-17 18-23 224 Totals 
Speed (kt) > Spells Hours 
Summer 25 515 33 25 8 5 4 882 2992 
30 168 C 12 5 I 2 I 997 
35 63 : 5 I 260 
40 22 I 84 
45 10 27 
50 4 12 
55 I 4 
Winter f 9 p 10 O19 
30 733 6 4 430 
287 1 418 


421 
30 96 
12 9 27 
6 9 
3 + 


gol 
350 
40 
16 


| 
3 
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TABLE III—TOTAL NUMBERS OF HOURS WITH MEAN WIND SPEEDS EQUAL TO OR 
GREATER THAN STATED VALUES IN EACH MONTH DURING THE 14 YEARS 1957-70 AT 
LERWICK, IRRESPECTIVE OF DIRECTION 


Mean speed equal to or greater than 
25 kt gokt 35kt gokt 45kt 50kt 55 kt 
Numbers of hours 

January 1 981 919 317. ~=«109 26 
February 464 642 201 89 29 

March 998 949 210 40 8 

April 799 «6295 8985 7 

May 445 173 46 10 

June 584 154 7 

July 259 41 oe 

August 265 71 20 


September 640 263 87 20 
October 309 522 86184 22 
November 302 422 120 2 
December 965 976 =386 9 


Year 13 O11 5427 1678 505 123 39 13 4 


have been very effective in generating high waves either because the higher 
speeds were of limited duration or because there were marked changes in wind 
direction. An alternative selection of notable storms could be made on the basis 
of a high average speed over a longer period, say 40 kt or more over 12 hours, 
combined with a limited direction range, say 40 degrees or less. There were 13 
storms which met these criteria, seven of which already appear in Figure 1 
(numbers 2, 3, 6, 8, 10, 11 and 12) and the other six of which are shown in 
Figure 2. In both figures a horizontal line drawn beneath each relevant speed 
profile shows the 12-hour period giving the highest average speed, and is 
labelled to show the average speed and direction range concerned. 

Like the 15 storms of Figure 1, these 13 12-hour storms have no common 
feature, although the majority of the speed profiles show an increase followed by 
a decrease, and the majority of the direction profiles show a more or less steady 
veer with time. Thus it cannot be said that there is a typical storm profile. 


Extreme wind speeds. Table IV shows the highest mean wind speeds, 
reduced to the standard height of 10 m, recorded in each of the years 1957 to 
1970 for durations ranging from 1 hour to 72 hours and for direction ranges of 
40 degrees or less. In any one year the highest speeds are given for each duration 
irrespective of the storm in which they occurred, but as the times and dates of 
commencement are given it can easily be seen which spells occurred in the same 
storm. Out of the 14 years there were 9 in which hourly mean speeds of 50 kt or 
more occurred, the highest speed being 63 kt in January 1961. Nine years also 
produced 6-hour mean speeds of 4.5 kt or more, seven years gave 12-hour means 
of 40 kt or more while 11 years gave 24-hour means of 35 kt or more. Even over 
as long a period as 48 hours, eight of the 14 years gave mean speeds of go kt or 
more. 

Tables V and VI give the annual extreme speeds over 12-hour and 24-hour 
periods respectively, together with the highest speeds averaged over shorter 
periods in the same storms. As in Table IV the speeds were reduced to the 
standard height of 10 m and are for direction ranges of 40 degrees or less. 
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TABLE V——HIGHEST MEAN WIND SPEED IN KNOTS OVER A I2-HOUR PERIOD DURING 

WHICH THE DIRECTION RANGE DID NOT EXCEED 40° IN EACH OF THE YEARS 1957 TO 

1970 AT LERWICK, TOGETHER WITH MAXIMUM SPEEDS OVER PERIODS OF I HOUR 

AND 3, 6, AND QO HOURS IN THE SAME STORMS; ALL SPEEDS ARE REDUCED TO VALUES 
APPROPRIATE TO A HEIGHT OF IO METRES 


Duration (hours) 


6 9 12 Times of 
knots commencement 
49°2 44°0 1th 5 Feb. 
40°3 39°38 39°4 th 29 Dec. 
49°1 484 46:5 20h 27 Oct. 
40°6 38-9 oh 5 Apr. 


49't 463 42°9 7h 27 Jan. 
422 40°83 393 1gh 15 Feb. 
54°5 506 468 14h 26 Sept. 
40°7 40°2 39°7 12h 14 Mar. 
52°2 46°8 14h 28 Oct. 


47°0 448 1th 6 Sept. 
40°7 39°5 4h 3 Dec. 
46-2 44°7 8h 4 Feb. 
45°5 : 39°! 2th 28 Sept. 
40°5 37°8 2th 23 Apr. 


TABLE VI—HIGHEST MEAN WIND SPEED IN KNOTS OVER A 24-HOUR PERIOD DURING 

WHICH THE DIRECTION RANGE DID NOT EXCEED 40° IN EACH OF THE YEARS 1957 TO 

1970 AT LERWICK, TOGETHER WITH MAXIMUM SPEEDS OVER PERIODS OF I HOUR 

AND 3, 6, 9, 12, 15 AND 18 HOURS IN THE SAME STORM; ALL SPEEDS ARE REDUCED 
TO VALUES APPROPRIATE TO A HEIGHT OF 10 METRES 


Duration (hours) 


3 9 12 15 18 24 Times of 
knots commencement 
1957 39°4 6 385 381 37:0 36:3 36:7 20h 7Jan. 
1958 40°6 398 394 385 37°55 35°9 2th 28 Dec. 
1959 49°3 ‘1 48:4 465 845°6 38448 43:4 gh 27 Oct. 
1960 42°4 39°5 389 382 37:9 37° 13h 4 Apr. 


1961 ‘ 55°8 46°3 42°99 40°3 39°83 4I1°0 19h 26 Jan. 
1962 43°3. 42:2 408 393 37°99 36:1 33°5 toh 15 Feb. 
1963 45°09 43°09 41°2 398 38:3 37:2 36°6 15h 24 Dec. 
1964 413. 40°7 402 39°7) 39°2 38:7 37°7 gh 14 Mar. 
1965 443 41°5 406 388 37:5 36:3 32:8 18h 30 Oct. 


1966 46-3 45:0 43°83 41°8 39:9 38:2 36-0 15h_ 6 Sept. 
1967 45°0 42°55 409 39°! 37°7 37°4 35°6 5h 21 Mar. 
1968 470 46:2 45°7 44°7 43°7 42:3 37°7 3h 4 Feb. 
1969 39°3 39°90 389 6386 0383) 37°99 3771 2h 14 Dec. 
1970 413. 40° «39 389378 «36:9 635°2 32-7. ~=Igh 23 Apr. 


Amongst other things the data in Tables IV, V and VI may be used to esti- 
mate how extreme wind speeds fall off as the averaging period is increased be- 
yond one hour. Average speeds over the 14 years for each duration were com- 
puted from each table and expressed as ratios of the averages over one hour, and 
the results are given in Table VII. As was to be expected, the mean ratios 
derived from the same storms for each year, that is to say from Tables V and VI, 
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are somewhat higher than those based on the highest speeds irrespective of the 
storms in which they occurred, as given in Table IV, particuarly for the longer 
averaging periods. It is suggested that the ratios in line (a) of Table VII should 
be used when, given an extreme hourly mean speed such as might be computed 
from anemograph records or interpolated from a map of once-in-50-year hourly 
mean speeds, it is desired to estimate an extreme for a longer period, say 24 hours, 
having the same probability, it being understood of course that the longer period 
is one in which wind direction remains within a 40-degree range. Thus if the 
once-in-50-year hourly mean speed was 65 kt then the estimated once-in-50-year 
12-hour and 24-hour means would be about 65 x 0-82 = 53 kt and 65 x 0-71 
= 46 kt respectively. The ratios on lines (b) and (c) of Table VII on the other 
hand might be used when, given an extreme for one averaging period, it was 
desired to estimate probable extremes for other, longer or shorter, averaging 
periods in the same storm. Thus, given a highest 24-hour mean of 45 kt ina 
storm with no great direction change the probable highest 3-hour mean in the 
same storm would be about 45 x 0-97/0-80 = 55 kt; or given a highest hourly 
mean of 60 kt the probable 12-hour extreme in the same storm, assuming no 
great direction change, would be about 60 x 0-86 = 52 kt. Also shown on line 
(d) of Table VII are ratios previously estimated by Shellard and published by 
the Department of Energy,’ which are in very good agreement with those now 
derived from 24-hour storms at Lerwick. 


TABLE VII—RATIOS OF MAXIMUM SPEEDS AVERAGED OVER VARYING NUMBERS OF 
HouRS (V;) TO MAXIMUM SPEEDS OVER ONE HOUR (V,) AT LERWICK 


Period ¢ (hours) 


3 12 18 24 36 48 72 
Ratio ;/V, 
0°96 0°90 0-82 0°74 o-71 . 0°58 Or51 
0°97 0°93 0-86 
0°97 0°93 0°88 0-83 0°80 
0-96 0°93 0°87 0°80 


derived from annual extremes for each duration irrespective of the storm in which they 
occurred; 

from maximum 12-hour storms; 

from maximum 24-hour storms; 

for comparison—values as given by Shellard in reference 1. 


The annual extreme wind speeds presented in Tables IV, V and VI may be 
fitted by extreme-value distributions of the Gumbel type, so providing estimates 
of extreme wind speeds for different averaging periods (durations) likely to be 
exceeded on average only once in, say, 50 or 100 years. Such estimates have 
been computed for average recurrence periods of 10, 20, 50, 100 and 200 years 
and are presented in Table VIII. Those in part (a) of the table were derived 
from Table IV and once-in-50-year extremes range from 71 kt for hourly means 
to 40 kt for 72-hour means. Those in parts (b) and (c) were derived from Tables 
V and VI and give for the same recurrence periods the probable extreme speeds 
in 12-hour and 24-hour storms respectively. In all cases direction ranges are 
assumed to be 40 degrees or less. 
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TABLE VIII—ESTIMATED EXTREME WIND SPEEDS AT LERWICK FOR VARIOUS 
DURATIONS AND AVERAGE RECURRENCE PERIODS 


(a) Estimated from data in Table IV, i.e. from highest speeds irrespective of 
the storms in which they occurred. 


Duration (hours) 
Average 


recurrence 

period 
I 3 6 12 18 24 36 48 72 
Extreme wind speed in knots for durations shown 


10 years 61 59 48 37 33 
20 years 65 63 50 40 36 
50 years 71 68 44 40 
100 years 75 72 47 43 
200 years 79 «676 49 46 


(b) Estimated from data in Table V, i.e. from highest speeds in 
maximum 12-hour storms. 


Duration (hours) 
Average 


recurrence 
period 
I 3 6 9 12 
Extreme wind speed in knots for durations shown 


10 years 61 57 48 
20 years 66 61 ; 50 
50 years 72 67 53 
100 years rai 71 56 
200 years 82 75 9 58 


(c) Estimated from data in Table VI, i.e. from highest speeds in 
maximum 24-hour storms. 


Duration (hours) 
Average 
recurrence 
period 


I 3 6 9 12 115 138 2% 
Extreme wind speed in knots for durations shown 


10 years 52 47 45 44 42 
20 years 55 49 47 45 
50 years 59 52 49 48 
100 years 62 54 51 49 
200 years ‘ 65 5 56 53 «+5! 


Application to wave prediction. As mentioned earlier, the original com- 
puter tabulations gave for each 30-degree direction range the numbers of spells 
with speeds equal to or greater than 20, 25, 30 etc. kt, and with durations of 
I, 2, 3, .- . 35 hours and 36 hours or more. By employing the wave-prediction 
technique of Darbyshire and Draper,? modified by Draper,’ these frequencies 
may be converted into wave-height frequencies. The relationship is not a 
simple one since wave height depends on both wind speed and duration, but the 
procedure used is illustrated in Table IX. In this table the numbers of 3-hour 
spells with speeds exceeding the values in column 1 and with directions in the 
range 230-250° is given in column 2. In column 3 the speeds of column 1 were 
adjusted to represent more closely the probable mean speeds over the open sea 
rather than those measured at Lerwick itself. The correction applied was a 
straightforward increase of 10 per cent, which figure was arrived at as follows. 
At coastal anemograph stations and in strong winds the average ratio of the 
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maximum gust speed to the maximum hourly mean speed, G (3-s, 60-min) has 
been shown by Shellard‘ to be about 1-5. Measurements over the sea by 
Goptarev,® Dorrestein® and Walden’ suggest that G (3-s, 10-min) is no more 
than 1-3 and this corresponds to a value of G (3-s, 60-min) of about 1°37. Also 
recent measurements on fixed gas-production platforms in the southern North 
Sea have shown G (3-s, 60-min) to be about 1-2 for anemographs at a height of 
80 m or so. When speeds are reduced to the standard height of 10 m using 
appropriate power-law formulae (exponent 0-12 for hourly means and 0-06 for 
gusts) this too gives a ratio of about 1-37. On the assumption that in strong 
winds maximum gust speeds will be much the same over the open sea as they are 
on nearby coasts, both being mainly dependent on the gradient wind speed, 
maximum hourly mean speeds over the sea will be about 1-50/1-37 or about 1-10 
times those on nearby coasts. 

Next it was necessary to decide which set of wave-prediction graphs provided 
by Darbyshire and Draper should be used, those for oceanic waters or those for 
coastal waters. On the advice of one of the authors (L. Draper, personal com- 
munication) the oceanic-waters graphs were used as being more appropriate to 
the Shetland area. It should be mentioned that the oceanic-waters graphs were 
derived from measurements made at the ocean weather stations ‘I’ and ‘J’, that 
is to say from wind speeds measured over the open sea. Had the coastal-waters 
graphs been used it would have been more correct to use the wind speeds as 
measured at Lerwick because these graphs were based on wind measurements 
made at coastal stations. The wave-prediction graphs provide estimates of (a) 
maximum wave height in feet during a 10-minute period and (b) significant 
wave period in seconds, for various combinations of wind speed and duration (or 
fetch). The appropriate values of these items are given in columns 4 and 5 of 
Table IX, wave heights being converted to metres. Columns 6 and 7 give the 
numbers of waves in the storm (in this case in three hours) and in the 10-minute 
period respectively, and column 8 gives the corresponding wave-height factors 
F, and F,, obtained from a diagram given by Draper (Figure 2 of reference 3) 
and their ratio. This diagram gives the ratio of maximum wave height to root- 
mean-square wave height for various numbers of waves. Since this ratio 
increases with number of waves, the ratio F,/F, represents the greater chance of 
a number of component waves getting into phase during the whole storm than 
in the 10-minute recording period on which the wave-prediction graphs are 
based ; thus the last column of Table IX gives values of Hmax (10-min) X F,/F,, 
the estimated maximum wave height in the storm. 

Hence, from each set of wind-speed frequencies, one for each direction range 
and duration, a set of highest-wave-in-storm frequencies can be obtained. 
These frequencies, derived from storms of various durations, can then be com- 
bined to give overall frequencies of waves exceeding various heights from each 
direction. This was done by plotting each set of frequencies against wave 
heights, for example column 2 against column g of Table IX, on log-linear 
graph paper and then reading off, by interpolation, the frequencies of wave 
heights at intervals of 14 metres from 4} metres upwards. These were then 
summed to give the results presented in Table X. 

In the process of deriving Table X some information was obtained which 
related highest storm waves to storm duration. This is summarized in Table XI 
which gives numbers of storms in which the predicted highest wave exceeded 
stated heights, arranged according to storm duration, all wind directions being 
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combined. It will be seen that all but one of the 18 storms in which the predicted 
highest wave was 15m or more in height had durations of less than 10 hours, that 
all five storms in which the predicted highest wave was 18 m or more in height 
had durations of less than 1o hours and that all three storms in which the pre- 
dicted highest wave was 21 m or more in height had durations of between 2 and 
6 hours. It appears that in the Shetland area the highest waves tend to be 
associated with rather severe gales of relatively short duration. Over the 14-year 
period there were only two storms having a duration of 24 hours or more and 
giving predicted highest waves of 10} m or more and only one storm of duration 
36 hours or more giving a predicted highest wave of 9 m or more. 

It should be noted that the predicted wave heights given in Tables X and XI 
take no account of swell. However, according to an analysis of wave data for 
ocean weather station ‘I’ (59°N, 19°W) by Hogben® the mean underlying swell 
in that area of the North Atlantic has a height of about 2 metres. The addition 
of an average swell wave of 2 m to wind waves of 6 m, 12 m and 24 m would give 
resultant waves of only 6-3, 12:2 and 24-1 m respectively. Even an exceptionally 
heavy swell of, say, 6 m when combined with sea waves of 12 and 24 m would 
give resultant waves of only 13-4 and 24-7 m respectively. 

Returning to Table X, it should be stressed that too much significance should 
not be attached to the actual frequencies given there, or in Table I, bearing in 
mind how they were obtained. Each 30-degree range was treated separately so 
that if, in a given storm, the wind direction veered or backed into an adjacent 
sector, that particular storm would have been terminated, even though the 
overall direction range may have been small enough for the storm to have been 
quite effective in developing waves. Also, the greater the duration of a storm the 
more likely would it be to be terminated in this way, thus being counted as two, 
or more, storms of lesser duration. However, the relative significance of the 
frequencies as given will not have been much affected by these occasional 
happenings and it is thought that the relative frequencies of waves of different 
heights from various directions are adequately represented by Table X. 

The table shows that the three highest waves, those of 21 m or more in height, 
came from directions between 230 and 280°, whereas those of 15 m or more, 
averaging just over one per year, came either from between 170 and 280° or 
from 350-010°. Between the directions 050° and 130° waves of 12 m or more in 
height appear to be unlikely in the Shetland area. 


Probable extreme wave heights. Finally, the Darbyshire/Draper wave- 
prediction technique may be applied to the extreme-wind-speed estimates given 
in Table VIII to provide estimates of the highest waves to be expected on aver- 
age only once in 10, 20, 50, 100 or 200 years. This has been done for an average 
recurrence period of 50 years using the appropriate wind speeds taken from part 
(a) of Table VIII and the results are given in Table XII. The values of Hmax 
(10-min) and of wave period enclosed in brackets are values whose derivation 
necessitated some limited extrapolation of the Darbyshire and Draper graphs. 
Since some preliminary calculations indicated that the predicted extreme wave 
heights increased rather quickly to a maximum value as the storm duration in- 
creased, interpolated extreme wind speeds for durations of 2, 4 and 5 hours were 
also used and the results are included in the table in view of the fact that the 
maximum value and the duration at which it occurs are of special interest. 
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The maximum once-in-50-year wave height of just over 35 m (about 115 ft) is 
associated with a once-in-50-year mean wind speed of 72 kt having a duration 
of about 4 hours. The corresponding significant wave period is about 19 
seconds. This maximum wave height compares with a once-in-50-year value 
of between 30 and 33 m to the west of the Shetlands taken from the map of 50- 
year design wave heights prepared by Draper and published by the Department 
of Energy.! Considering the various uncertainties involved this may be re- 
garded as a satisfactory degree of agreement. 


Acknowledgements. The author wishes to thank Mr R. Anderson-Jones 
who designed and wrote the necessary computer programs, and Mr L. Draper, 
Institute of Oceanographic Sciences, who kindly read a first draft of the paper 
and made some helpful suggestions. 


REFERENCES 


London, Department of Energy. Guidance on the design and construction of offshore 
installations, Annex to Section 2. London, HMSO, 1974. 

DARBYSHIRE, MOLLIE and DRAPER, L.; Forecasting wind-generated sea waves. Engineering, 
London, 195, 1963, pp. 482-484. 

DRAPER, L.; The analysis and presentation of wave data—a plea for uniformity. Pro- 
ceedings of the roth Conference on Coastal Engineering, Tokyo, 5-8 September 1966, 
Vol. 1. New York, American Society of Civil Engineers, 1967, pp. 1-11. 

. SHELLARD, H. C.; Extreme wind speed over the United Kingdom for periods ending 1959. 
Met Mag, London, 91, 1962 pp. 39-47. 

GOPTAREV, N. P.; Some results of wind profile measurements on marine oil drilling plat- 
forms. Trudy Gos Okeanogr Inst, Moscow, 36, 1957, pp. 128-202 (in Russian) as quoted 
by H. U. ROLL in ‘Physics of the marine atmosphere’, New York and London, Academic 
Press, 1965, pp. 168-170. 

). DORRESTEIN, R.; Unpublished internal report of the Koninklijk Nederlands Meteorologisch 
Instituut, De Bilt, the Netherlands, 1963. 

WALDEN, H.; Windmessungen im 10-sec-Abstand auf Forschungsschiff ‘“‘Meteor’’. Met 
Rundsch, Berlin, 24, 1971, pp. 27-29. 

HOGBEN, N.; Measured wave heights and wind speeds at weather station ‘India’ in the 
North Atlantic. Mar Obsr, London, 39, 1969, pp. 190-193. 


551.577-2(429) :551.589.1 


THE DISTRIBUTION OF RAINFALL OVER SUBCATCHMENTS OF 
THE RIVER DEE AS A FUNCTION OF SYNOPTIC TYPE 


By C. A. NICHOLASS and T. W. HARROLD 


Summary. Two years of data from a network of 63 tipping-bucket rain-gauges distributed 
over the 1000-km? catchment area of the upper portion of the River Dee have been analysed 
to derive relationships between the rainfall over subcatchments of area typically 60 km? (Rs) 
and that averaged over the entire area of the network (R). It is shown that the ratio Rs/R 
is dependent on the synoptic type and surface wind direction. Using these two parameters as 
inputs, the distribution of storm rainfall within the upper portion of the Dee Catchment could 
be forecast with adequate accuracy from an accurate quantitative forecast of the rainfall over 
the whole catchment area. The extent to which this conclusion can be applied to other areas 
is discussed. 


Introduction. A network of 63 tipping-bucket rain-gauges is operated by 
the Water Data Unit and the Welsh National Water Development Authority 
as part of the Dee Weather Radar Project (see for example Harrold, English and 
Nicholass'). These gauges provide data in 15-minute time periods over a 1000- 
km? catchment area of the upper portion of the River Dee in North Wales. 





Meteorological Magazine, 104, 1975 209 


Two years of the data obtained have been used to determine relationships be- 
tween subcatchment and catchment areal rainfalls for different synoptic 
weather types. The subcatchments vary in size from 20 km? to 104 km? (see 
Figure 1). 

Such relationships may be of use to the forecaster and hydrologist, since they 
would enable forecasts of rainfall amounts over large areas, such as those which 
will be available eventually from computer models, to be used in estimating 
rainfall over the much smaller areas which are of importance to hydrologists, 
particularly in the efficient operation of regulating reservoirs. 

Although the data demonstrate only the relationships which exist within the 
Dee Catchment, ways in which similar climatological rules could be investigated 
in other river catchment areas are discussed. 


Analysis and results. For the two-year period ending April 1974 the 
synoptic type and surface wind were classified by reference to the Daily Weather 
Report. The surface wind velocity over the Dee Catchment was estimated from 
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4°W — 1000 ft (305 m) 
Irish Sea 2000 ft (610 m) 




















20 km 
J 




















FIGURE I—MAP SHOWING THE LOCATION OF THE RIVER DEE CATCHMENT AREA 
AND THE 16 SUBCATCHMENTS REFERRED TO IN THE TEXT 
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nearby observations in four categories of wind direction, N to E, E to S, S to W, 
and W to N, and three categories of speed, 0-15 kt, 15-25 kt, and more than 
25 kt. (The possible situation that all of the relevant surface winds were exactly 
south or west for example, thus making classification difficult, did not occur in 
this analysis.) 

Early in the analysis it was found that when the surface wind speed was more 
than 15 kt the actual speed was not a significant factor in determining the 
distribution of rainfall within the Dee Catchment. (Not enough heavy rainfalls 
occurred when the wind speed was less than 15 kt to allow a decision on the 
importance of light wind speeds to be made.) Therefore only surface wind 
direction on a four-point compass (SW means winds between S and W) is 
discussed here. 

Six synoptic types were classified: pre-warm-front; warm sector (including 
rain ahead of cold fronts); post-cold-front; occlusions; cyclonic rain not 
associated with well-defined frontal systems; and showers (with or without 
thunder) not associated with any other type. Periods during which the synoptic 
type and surface wind direction were constant are referred to as weather types 
in the following; they lasted between 3 and 36 hours approximately. For each 
weather type the areal rainfall for the whole of the Dee Catchment and for each 
of 16* subcatchments was computed by the method described by English.? 
Occasions when snow may have fallen or may have been lying in the gauges 
have not been included in the analysis. 

Graphs were plotted of catchment rainfall against subcatchment rainfall for 
the different weather types. They showed that a linear relationship existed 
between the two parameters, so a line of best fit, correlation coefficient and 
standard error of the estimate were calculated for each subcatchment and 
weather type. 

Examples of the data from two subcatchments for warm-sector rain with 
south-west winds are shown in Figure 2. The total number of rain periods in 
this category was 52. These graphs and those for each subcatchment for this 
weather type are summarized in Table I. The scatter about the line of best fit is 
expressed in two ways. In column 4 the standard error of estimate of Rs is shown 
in millimetres. In column 5 the error is expressed as a percentage of the mean fall 
over each subcatchment calculated from the 52 occasions. Averages of the 
errors shown in columns 4 and 5 are 2:3 mm and 33 per cent respectively. 
These values are indicative of the errors to be expected in a forecast of sub- 
catchment rainfall, and they show that in the particular type of situation under 
consideration, the distribution of the rainfall over the subcatchments could be 
predicted with a fairly high degree of accuracy, provided that the rainfall over 
the entire catchment could be correctly forecast. 

Table II summarizes some of the statistical parameters for the 10 weather 
types in which sufficient rain fell for significant regression equations to be 
calculated. Data for three subcatchments, which are representative of the 15, 
are presented. These results show that the scatter in the relationships in frontal 
rains is generally similar to those in Table I. Not surprisingly the scatter is 





* Note: The results from subcatchment 12 will not be discussed since the small number of 
gauges in the area made it difficult to obtain accurate areal estimates over the subcatchment; 
the catchment estimate should not be significantly influenced by the sparser network in 
this region. 
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TABLE I—WARM-SECTOR RAINS WITH SOUTH-WEST WINDS 


Standard error of 
estimate of Rs a 
Slope of line of Correlation (Ss) Ss/Rs 
Subcatchment best fit Rs/R coefficient millimetres per cent 

I 0°66 0°85 2°3 
0°73 0°88 2°2 
I-O1 0°95 1-8 
1°35 0°97 2°0 


OO ND 


1°76 0:89 52 
1°74 0°92 4°2 
1°55 0-96 2°5 
1°44 0°97 2°2 


Cnr OU 


1°12 0°99 
1°04 0°96 
o-71 0°88 
0°52 0°86 


-N- © 


wd BUGS 


0°59 0-89 
0°77 0-96 
0°77 0°95 


~~ 


Average values 
* Subcatchment 12 has been omitted. 


Ls) 
=) 


largest in the eight occasions of showery rain which did not fit into one of the 
other classifications. In this class there appeared to be little or no systematic 
effect. 

The results also show that the orographic effects are strongly dependent on 
weather type. The variability of rainfall on the scale of subcatchments is most 
marked within warm sectors and least marked ahead of warm fronts when the 
surface wind is south-westerly (see Table II). These findings are consistent with 
the conclusions of Browning e¢ alii,? which were based on a very limited, but 
extensively analysed, number of cases. 

The preceding statistical results have been computed for rainfall totals from 
periods of constant weather type. Falls from a few of the wetter storms shown in 
Figure 2 have been subdivided into hourly totals to investigate the extent to 
which these climatological rules are applicable to shorter periods within the 
storms. An example of the results obtained is shown in Figure 3. The line of 
best fit has a slope of 1-57 with a standard error of estimate(S) of 0-82 mm and an 
S/Rs of 30 per cent. So even on this short time-scale the distribution of rainfall 
can be predicted to a quite high degree of accuracy, provided that the large- 
scale rainfall can be forecast perfectly over the whole catchment. 


Implications. The results of this investigation show that in general the 
amount of rain falling in a given synoptic situation over any subcatchment of the 
upper portions of the River Dee is closely related to the rainfall over the entire 
catchment. That is to say, given the large-scale (synoptic-scale) precipitation, 
the topography is by far the most important factor in controlling the distribution 
of rainfall in this area. The finding has important implications in forecasting 
rainfall amounts over subcatchments for hydrological purposes. Evidently the 
accuracy of forecasts of rainfall on the scale of subcatchments depends primarily 
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FIGURE 2—RAINFALL IN AREA 6 (4) AND AREA 9 (©) PLOTTED AGAINST AVERAGE 
RAINFALL FOR ENTIRE AREA OF THE NETWORK (#) FOR WARM-SECTOR RAINS 
(AND RAIN AHEAD OF COLD FRONTS) WITH SOUTH-WEST WINDS 
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FIGURE 3—HOURLY RAINFALL IN AREA 6 (R,) PLOTTED AGAINST HOURLY 
AVERAGE RAINFALL (R) FOR A SAMPLE OF THE WARM-SECTOR RAINS 
WITH SOUTH-WEST WINDS 


on the accuracy with which the rainfall amounts over the entire catchment of 
10° km? can be forecast. With this proviso that the large-scale rainfall can be 
accurately forecast, the distribution of rainfall within the catchment is reason- 
ably predictable, particularly in frontal situations which produce by far the 
largest portion of the total rain in the area. The main exception occurs in 
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showery conditions. In addition, a mesoscale rain band not related to the topo- 
graphy (at least in its immediate vicinity) may occasionally dominate the frontal 
rainfall pattern for a period. Further research is needed in order to identify in 
advance occasions when such exceptions might occur. 

The data were gathered over one particular hilly area. It is not known to 
what extent analogous synoptic climatological relationships apply in other hilly 
areas. However, the topography of the catchment of the Upper Dee is quite 
complex, with mountain ranges at several different orientations, so that the 
effects of the topography are complicated. Hence there is every reason to 
expect that similar well-defined relationships exist in other hilly areas. Un- 
fortunately the data needed to determine such relationships cannot easily be 
obtained, since rainfall totals over small areas are required for periods of 
‘constant’ synoptic type. Ways in which such data might be obtained are: 


(a) to use existing ; raphic rain-cauge data. However, the density of 
autographic gauges is t arse in almost 7 hilly regions to determine sub- 


catchment rai 


(b) to obtain additional rain-gauge data. This could be done by greatly 
increasing the density of autographic rain-gauges in areas of interest. However, 
the cost of collecting and analysing the data from an adequate network of such 
gauges in hilly terrain is high—f40 000 per annum for the Dee network. An 
alternative, simpler, means of obtaining the necessary data is to use quantitative 
rainfall measurements derived from weather radar (Harrold et alii). The 
feasibility of doing this will be investigated using a mini-network of three 
quantitative weather radars during 1975 (Taylor and Browning‘), but this net- 
work will only be operated occasionally for special research studies. A routine 
operational network of radars, such as discussed by Dee Weather Radar Project® 
would be required if radar were to provide the amount of data needed to deter- 
mine climatological relationships over an extensive region; and 


(c) to use a numerical model of the topographic effects on rainfall. Such a 
model has been described by Collier,* and Table III shows ratios of R,/R over 
subcatchments of the Dee derived from this model in moving baroclinic dis- 
turbances, together with those from the rain-gauge network, for south-westerly 
winds. The average difference in R,/R between these entirely different tech- 
niques is 21 per cent. However, it is not yet certain to what extent the model can 
be applied in other hilly areas; some independent data are required to investi- 
gate this. 

It must be stressed that synoptic climatological relationships between rainfall 
amounts over subcatchments (typical area 60 km?) to that over 1000 km? will 
only be of practical value if the rainfall amount over the 1000 km? can be 
accurately forecast. Possible methods of achieving reliable forecasts on the scale 
of 1000 km? are: 


(a) to use existing numerical techniques with a smaller (mesoscale) grid 
length. It is not yet known what the smallest scale of accurate forecasts is 
but it seems improbable that it will be less than 10% km? in hilly terrain; 


(b) to use climatological relationships similar to those described on page 210 
to link the smallest scale accurately forecast in (a) to the scale of 10% km? 
(and hence 10? km?). Suitable data do not exist at present for the determination 
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of relationships on this larger scale, but some should be forthcoming from the 
mini-network of research radars; and 


(c) to use a parameterization such as Collier’s* model to estimate the precipi- 
tation on a scale of 103 km? from an input of the larger-scale wind and humidity 
field, derived for instance from (a). This has been shown to be a successful 
technique over the Dee Catchment in moving baroclinic systems, using actual, 
rather than forecast, input parameters. 


TABLE III—COMPARISON OF PREDICTED AND ACTUAL CLIMATOLOGICAL RATIOS FOR 
WARM-SECTOR RAINS WITH SOUTH-WEST WINDS 


R;/R computed by Rx/R measured Computed ratio 
Subcatchment Coollier’s* model by rain-gauges = measured ratio 
per cent 


0°83 0:66 126 
0°75 0°73 103 
0°92 1-O1 g! 
1°41 1°35 105 
1°83 1°76 104 
167 1°74 96 
1°50 1°55 97 
1-08 1°44 75 
0°92 1°12 82 
125 1°04 120 
1-08 0-71 152 
067 0°52 129 


0°83 0°59 141 


0°92 0°77 11g 
1°16 0°77 I5I 


Mean absolute percentage difference = 21% 


Conclusions. The analysis has shown that in most circumstances over the 
upper portion of the River Dee forecasts of the rain over subcatchments of area 
20-100 km? could be forecast if there were a means of forecasting accurately over 
an area of 1000 km?. It is considered that a similar conclusion would apply 
to other hilly terrain. Thus, if a method of forecasting rainfall over the larger 
area were developed this would also enable the hydrological requirement of 
quantitative rainfall forecasts over the smaller areas to be met. 

However, there are several other forecasting requirements, for instance 
variations of rainfall intensity with time and also localized storms and stationary 
mesoscale rain bands which are not handled by the climatological rules. To 
meet these full requirements, it will probably be necessary to use results of 
the type presented in this paper in conjunction with real-time data from an 
operational network of weather radars. 
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processed by the Welsh National Water Development Authority and the Water 
Data Unit, Reading as part of the Dee Weather Radar Project. The authors 
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CONFERENCE ON ‘ENGINEERING HYDROLOGY TODAY’, 
LONDON, 18-20 FEBRUARY 1975 


By R. MURRAY 


A conference on ‘Engineering Hydrology Today’ was held at the Institution of 
Civil Engineers in London from 18 to 20 February 1975 under the joint sponsor- 


ship of the Institution of Civil Engineers, the Institute of Hydrology, and the 
organizing committee of the International Hydrological Decade. The confer- 
ence set out to review the contributions to the International Hydrological 
Decade (1965-74) made by the United Kingdom, with particular emphasis 
on results of relevance to engineering. 

The conference was opened at 2 p.m. on the 18th by Lord Nugent, chairman 
of the National Water Council, who spoke briefly about the reorganization of 
the water industry which took place on 1 April 1974 as a consequence of the 
Water Act 1973. It was appropriate that Sir Norman Rowntree, erstwhile 
Director of the former Water Resources Board, which ceased to exist on 31 
March 1974 as a result of the Act, should give the last talk on ‘Summing-up and 
a look to the future’. Between the contributions of Lord Nugent and Sir 
Norman Rowntree, the conference was organized into six sessions which dealt 
with Organization, Instrumentation and Techniques, Meteorology, Flow 
Models, Flow Frequency Estimation, and Storage. There were three papers 
directly concerned with meteorology, namely (a) ‘Determining precipitation, 
evaporation and soil moisture’ by Dr. J. C. Rodda (Water Data Unit) and 
Mr J. F. Keers (Meteorological Office), (b) “The variability of precipitation 
and evaporation’ by Mr J. F. Keers and Dr J. C. Rodda and (c) ‘Estimation 
of irrigation needs’ by Mr B. G. Wales-Smith (Meteorological Office). All three 
papers were well received by an audience of nearly 200. 

Dr Rodda, who introduced the first paper, demonstrated progress by the 
developments in measuring rainfall by radar (a major part having been played 
by the Meteorological Office team in the co-operative Dee Weather Radar 
Project) and in neutron probes, although Dr Penman was quick to point out 
that neutron probes were not novel to the Decade. Dr Rodda had to confess, 
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however, that rainfall continued to be measured for the most part by conven- 
tional gauges, but mentioned the Kew gravimetric gauge as a standard against 
which other types of gauge could be compared. There was lively discussion 
about the accuracy with which rainfall was measurable in different regions: 
in particular Mr Reynolds of the North of Scotland Hydro-electric Board 
thought that Dr Rodda was much too pessimistic in quoting a 20 per cent 
difference between readings in hilly terrain from the Meteorological Office 
Mk 2 rain-gauge sited in the standard way and the pit gauge with its rim at 
general ground level. Support for the accuracy of readings from the Mk 2 
gauge, properly sited, came from Dr Penman and others. Nevertheless, it was 
felt that the problem of accurately measuring rainfall was still with us, especially 
over difficult terrain, although it was recognized that useful areal estimates 
could be obtained by radar. Snow measurement continued to be an unsolved 
problem. The need to improve the reliability of recording gauges, which 
generally have too high a failure rate, especially in freezing weather, was 
stressed by several delegates. 

Mr Keers introduced the second paper but concentrated mainly on rainfall 
variability. A useful account was given of the different time and space scales 
of rainfall variability, with some reference to the precipitation mechanisms and 
underlying causes. He drew attention to some results of the work done by 
British radar meteorologists, notably Browning and Harrold, during the Decade, 
and referred briefly to the Meteorological Office’s ‘Project Scillonia’ and the 
recent Global Atmospheric Research Programme Atlantic Tropical Experiment 
(GATE) during the summer of 1974, in which the United Kingdom, notably the 
Meteorological Office through the Meteorological Research Flight, played an 
important part. Mr Keers also touched upon the recently completed Flood 
Studies Report, but detailed discussion of this report was reserved for a conference 
at the Institution of Civil Engineers in May 1975. 

Mr Wales-Smith surveyed a very wide field in his written paper but sensibly 
curtailed his spoken presentation and concentrated for the most part on studies 
in hand within the Meteorological Office. These studies concerned improve- 
ment of the ‘Estimated Soil Moisture Deficit? (SMD) bulletin by means of the 
incorporation of a more realistic land-use model and by an improv ‘eo computer- 
based system for accessing all types of up-to-date data needed for full exploita- 
tion of the Penman formula which constituted the scientific basis of the Meteoro- 
logical Office’s SMD bulletin. Briefreference was also made to rainfall deficiency 
studies and to other in ns. All such work should make it feasible to 
monitor soil moistu its mo tly thr hout the country and 
should help therefore in the monitoring of irrigation needs, especially if weather 
forecasts were sensibly 1 ine with the latest SMD information. 

7 e tw aner: ; i a} 7 1 \7 ’ “4 were 1°. ~ussed together 
and | ivel ind engineers. At one 
point there was danger of the discussion st ing too far into the field reserved 
for discussion at the « >on the Flood Studies Report in May 1975, but this 
clearly indicated the need felt by engineers for inf rmation on flood- -producing 
rainfall, especially rainfall of short duration, information which is important in 
planning drainage systems. The need for lnGeesiell advice on rainfall to 
engineers engaged on projects in hilly areas was mentioned. The lack of 
information on the variability of evaporation, comparable to the information on 
the variability of precipitation, was also mentioned in the discussion. Several 
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speakers clearly indicated the importance of work on evaporation and on soil 
moisture deficits to water resources management and indeed to civil engineers. 
Mr N. J. Cochrane (Sir William Halcrow and Partners) stressed the need 
for more information about evaporation and soil moisture conditions. Inci- 
dentally Mr Cochrane mentioned that neutron probes had been used—without 
much success owing to difficulties of maintenance and power supply—in two 
road-building projects in Africa many years ago. However, there was little 
doubt that under reliable management neutron probes were now reliable 
instruments. 

Other important papers by hydrologists and engineers on ‘Catchment 
modelling to estimate flows’, ‘Open-channel hydraulics’, ‘Flow frequency 
estimates’, ‘Groundwater yield estimates from models’ and ‘Assessment of 
surface water sources’ were not so directly of interest to meteorologists, although 
the meteorological input (precipitation and evaporation) is of great importance. 
It was interesting to a meteorologist to learn that some of the engineers and 
practical hydrologists were rather sceptical about the value of complicated 
mathematical models and seemed to prefer empirical models, based on data and 
experience. There was also a feeling of distrust among some delegates of synthetic- 
data generation, especially in view of the amount of rainfall data, and to a 
less extent stream-flow data, available in the United Kingdom. 

The conference was particularly appreciated by the three Meteorological 
Office representatives who were able to meet engineers and hydrologists and to 
learn about some of their practical problems and their meteorological require- 
ments, as well as to enjoy accounts of some British achievements in hydrology 
during the recent Decade. 


PUBLICATIONS RECEIVED 
The following have been received from the Meteorological Institute of the 
University of Thessaloniki: 
Meteorologika 36: Nature of the diurnal variation of atmospheric pressure in Thessalonikt. 
By T. J. Makroyannis. 1974. 
Meteorologika 37. Cooling power and weather types in Thessaloniki. By Chr. J. 
Balafoutis. 1974. 


Meteorologika 38: On the effect of ground relief upon sunshine duration on Mount 
Olympus—Greece. By G. C. Livadas and V. A. Semertzidis. 1974. 


Meteorologika 39: On the annual variation of air temperature in Thessaloniki. By A. A. 
Flocas and A. Arseni-Papadimitriou. 1974. 


Meteorologika 40: Contribution to the study of air temperature in Thessaloniki. By 
A. Arseni-Papadimitrou. 1974. 


Meteorological Magazine: price increases 


As from July 1975 the price of an issue of the Meteorological Magazine will 
be 4op and the annual subscription will be £5.46 including postage. 
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NOTES AND NEWS 


Retirement of Mr J. M. Craddock 


On 4 June 1975 Mr James Marston Craddock retired from the Meteorological 
Office after 33 years’ service. For the previous 16 years he had been a Special 
Merit Senior Principal Scientific Officer concerned with the application of 
statistics and computers to long-range weather forecasting. 

Mr Craddock joined the Office in 1942 on secondment from the Inland 
Revenue Department, having previously taken a first-class degree in mathematics 
at Cambridge University where he was a scholar of Magdalene College, and 
spent the next five years of his career as a forecaster. After a short spell at 
Prestwick he was posted to the upper-air bench at Dunstable. He was then 
commissioned as a Flight Lieutenant in the Royal Air Force and posted to the 
Far East where he spent about a year in Singapore before moving to Butterworth. 

On demobilization in 1947, Mr Craddock, having decided that he would 
get more satisfaction from meteorological research than from returning to the 
Inland Revenue, was posted as a Senior Scientific Officer to Dunstable and was 
one of the first members of Dr Sutcliffe’s forecasting-research group in the 
Napier Shaw Building. Mr Craddock had clearly found his niche in life and, 
apart from a two-year spell in 1953-54 when he was posted to the Central 
Forecasting Office as a senior forecaster on his promotion to Principal Scientific 
Officer, he continued to work in the research branch connected with long-range 
weather prediction until his retirement. In 1959 he obtained a well-deserved 
special merit promotion to Senior Principal Scientific Officer. In addition he 
received the L. G. Groves and Darton Prizes in recognition of his scientific 
ability. 

Two of Mr Craddock’s major achievements in the Office have been the 
development of sound statistical methods for use in meteorological research 
and the application of computers to the operational and research work of the 
Synoptic Climatology Branch. In the statistical field Mr Craddock has made 
important contributions on topics such as the analysis of time series and the 
application of principal component analysis to meteorological problems. 
In the computer field he has been largely responsible for building up a large 
long-range data bank and for developing the METOCODE language which 
has considerably reduced the amount of programming effort required to permit 
statistical programs to be run on the computer. In addition, Mr Craddock 
has been a World Meteorological Organization consultant concerned with the 
collection, storage and cataloguing of meteorological literature and data. 

We all wish Mr and Mrs Craddock many years of happy retirement. 


F. H. BUSHBY 
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